Estimation of fluctuating magnetic flelds by an atomic magnetometer 
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We present a theoretical analysis of the ability of atomic magnetometers to estimate a fluctuating 
magnetic field. Our analysis makes use of a Gaussian state description of the atoms and the probing 
field, and it presents the estimator of the field and a measure of its uncertainty which coincides 
in the appropriate limit with the achievements for a static field. We show by simulations that the 
estimator for the current value of the field systematically lags behind the actual value of the field, 
and we suggest a more complete theory, where measurement results at any time are used to update 
and improve both the estimate of the current value and the estimate of past values of the B-field. 
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I. INTRODUCTION 

For decades superconducting quantum interference de- 
vices, SQUIDs, have been unrivalled as the most sen- 
sitive detectors of weak magnetic fields, and they have 
been applied in diverse scientific studies including NMR 
signal detection IM, visualization of human brain activ- 
ity and gravitational wave detection [3'j. Optically 
pumped atomic gasses offer an alternative means to de- 
tect weak magnetic fields via the induced Larmor preces- 
sion of the polarized spin component, and the possibil- 
ity to avoid the need for cryogenic cooling and the rela- 
tively high price of the SQUID devices, has spurred an 
interest in employing atomic magnetometers in medical 
diagnostics, see for example work on the mapping of hu- 
man cardiomagnetic fields, Combined with the fact, 
that the atomic based magnetometers may now reach su- 
perior field sensitivity Q, we are now highly motivated 
to investigate and identify the optimal performance and 
fundamental limits of such devices. 

Atoms constitute ideal probes for a number of physical 
phenomena, and since they are quantum systems, the sta- 
tistical analysis of measurement results has to take into 
account the very special role of measurements in quan- 
tum theory. In high precision metrology the aim is to 
reduce error bars as much as possible, and there is an ob- 
vious interest in making the tightest possible conclusion 
from measurement data. How to optimally prepare and 
interact with a physical system to obtain maximum infor- 
mation, and even the simpler task of identifying precisely 
the information available from a specific (noisy) detection 
record are not fully characterized at this moment. Large 
efforts are currently being made to combine techniques 
from classical control and parameter estimation theory 
with quantum filtering equations. 

In the present paper we analyze the case of atomic 
magnetometry with large atomic samples which are be- 
ing probed continuously in time with weak optical probes. 
Such systems allow a simplified treatment by means of 
Gaussian states and probability distributions, and we are 
able to give exact expressions for the probability distribu- 
tion of the probed magnetic field, i.e., our result is neither 



too weak nor too strong (no procedure exists by which 
further information can be extracted from the available 
data, and the actual value must agree with our estimate 
within the probability distribution). Our work is a gener- 
alization of previous work for static fields 0, 0| to the in- 
teresting case of time dependent fields. Time dependent 
fields were studied recently and we shall generalize 
that analysis and show that further improvements of the 
estimate of the field at a given time is possible by taking 
into account the detailed detection record at both earlier 
and later times. 

In Sec.m we describe the physical setup and the inter- 
actions between a single component B-field , the atomic 
system and the optical probe field. In Sec. Ill Bl we de- 
scribe the time evolution of the collective atomic quan- 
tum state during interaction with a given, parametrized, 
field B{t) and optical probing. In Sec. IIIII we treat the 
case, where the time dependent magnetic field is assumed 
to be a realization of an Ornstein-Uhlenbeck stochastic 
process, where the actual time dependence of B{t) is un- 
known by the experimentalist, and we provide the formal- 
ism for the optimum estimate for B{t) based on the noisy 
detection record for all earher times t' < t. In Sec. IIVI 
we show results of numerical simulations, suggesting that 
also the detection at later times t' > t improves our esti- 
mate of B(t), and we present ad hoc and precise analyses 
of the optimum estimator. Sec. concludes the paper. 



II. INTERACTION BETWEEN A LARGE 
ATOMIC SAMPLE, A MAGNETIC FIELD AND 
AN OPTICAL PROBE 

A. Physical System and Interactions 

We consider a gas with a macroscopic number TVat 
of atoms with two degenerate Zeeman states. Initially, 
all the atoms are prepared by optical pumping in the 
same internal quantum state, and all interactions are as- 
sumed to be invariant under permutations of the atoms. 
The dynamics is conveniently described by the collective 
effective spin operator J = f X]j with (Ti the Pauli 
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FIG. 1: A B-field component along the y-axis causes a Larmor 
precession of an atomic spin which is initially polarized along 
the a;-axis. The resulting z component of the spin is probed 
by the Faraday polarization rotation of a linearly polarized 
laser beam which propagates along the j/-axis. 



matrices describing the individual two-level atoms. The 
atoms are initially prepared by optical pumping such that 
their spin is polarized along the x axis, and we assume 
only a small depolarization during the interaction, so that 
the operator Jx can be well approximated by a constant 
number (Jx) — The other two projections of the 

collective spin, Jy and Jz, obey the commutation rela- 
tion [Jy,Jz] — ihJx and the resulting uncertainty rela- 
tion on Jy and Jx, i.e., on the number of atoms pop- 
ulating the (T^ and cr^ atomic eigenstates precisely re- 
flect the binomial distribution of atoms on these states. 
The commutator may be rewritten as [xatiPat] = * for 
the effective canonical position and momentum variables 
— — — — and the binomial population 



statistics of the collective states nicely maps to Gaussian 
probability distributions of x^t and Pat • 

We now imagine that the atomic sample is placed in 
a B-field directed along the y direction, which causes a 
Larmor rotation of the atomic spin towards the z axis. 
The setup is shown in Fig. During a time interval r, 
the z-component of the collective spin, represented by 
the operator Pat , evolves as 



Pat 



Pat - 



(1) 



where /i is give n by the magnetic moment /3, via = 
{l/h)(3yj {Jx) /h. By probing the value of Pat, we acquire 
information about the magnetic field. 

A number of papers have dealt with the interaction of 
atomic samples with an optical field, [6, 9, 10, 11, 12, ij 
Il3 |. A particularly convenient interaction occurs when 
the two atomic states interact with a linearly polarized 
non-resonant probe field. A linearly polarized field can 
be expanded on two circular components, and if they 
interact differently with the two internal atomic states, 
they experience different phase shifts depending on the 
atomic populations, which in turn translates into a (Fara- 
day) rotation of the polarization of the field. For use of 
the effective spin- 1/2 picture to the description of angu- 
lar momenta larger than 1/2, see for example We 
are interested in the case of a continuous beam of light in- 
teracting with the atoms, and this is effectively obtained 
by discretizing the field into a sequence of square pulse 
segments of light, see @. Since the light beam is linearly 
polarized along x, the Stokes operator for a short segment 
of the beam with A^ph photons has an x-component which 



TABLE I: Values of physical quantities used in the numerical 
examples of this paper. 



Physical quantity 



Value 



ABo = 1 pT 
A^at = 2 X 10^2 
$ = 5 X 10"s-i 
A = 2 mm^ 
A = 852 nm 
A = 10 GHz 



Initial uncertainty 
Total number of atoms 
Photon flux 
Intersection area 
Wavelength 
Detuning 

Atomic dipole moment d = 2.61 x 10~^^ Cm 

Coupling between atoms and B-field /i = 8.79 x 10'' 
Coupling between atoms and photons — 1.83 x 10^ s~^ 



behaves classically, {Sx) — ^^2^ , and the two remain- 
ing components fulfill the commutator relation of angu- 
lar momentum operators. Accordingly, for the scaled ef- 

s ^ 

we have 



fective variables Xph 



[xphiPph] = i and the initial coherent state of the field 
is a minimum uncertainty Gaussian state in these vari- 
ables. When a single light segment is sent through the 
atomic sample, the polarization rotation is quantified by 
the difference between the intensities of linear polariza- 
tion components at 45° and —45° angles with respect to 
the incident polarization, i.e., by the Sy (Xph) quantum 
variable. The field is not in an eigenstate of this quan- 
tity, and the measurement outcome is stochastic with a 
mean value governed by the mean rotation, i.e., the mean 
atomic population difference, and a variance given by a 
shot noise contribution from the photon number distri- 
bution, and by the variance of the atomic populations. In 
turn, the detection of the fight signal causes a back-action 
on the quantum state of the atoms, which, depending on 
the total number of photons detected, shifts the mean 
atomic populations and reduces their widths. 

As shown in R ef . [l4| and references therein, the polar- 
ization rotation described above is described by a simple 
effective Hamiltonian, and in the Heisenberg picture the 
atomic variables and the variables describing a single seg- 
ment of the light beam of duration r evolve as 



Xat a^at + Kv^Pph, 
Xph nVrPat + Xph, 



Pat Pat (2) 
Pph Pph (3) 



where the atom-light coupling strength k = 
AAca ^ ^a.tNp\y/T is 8, function of atomic parame- 
ters {d is the atomic dipole moment, lj is the photon 
frequency, A is the detuning of the light from atomic 
resonance, and A is the area of the the light field) and, 
notably, the total number of atoms A^at and the photon 
fiux $ = NpYi/r in the pulse segment. Our theory is to 
a large extent analytical, but in order to present some 
of the results in graphs, we shall adopt reaHstic values, 
summarized in Table 1, for a possible experiment with 
cesium atoms. 
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B. Gaussian state formalism 

In this section, we shall assume a magnetic field B{t) 
with an explicitly given time dependence. This analysis 
will be needed, when we proceed to simulate how an un- 
known field is estimated, since it is the actual realization 
of the field, that drives the atomic dynamics. The cur- 
rent section thus accounts for the "information available 
to the theorist", whereas the proceeding section will deal 
with the "information available to the experimentalist" 
having only access to the optical detection record. 

We place a spin-polarized atomic gas in the B-field, 
and probe it continuously by the polarized light beam. 
The Gaussian probability distribution (Wigner function) 
of the atomic variables and of the fight pulse prior to 
interaction, evolves into a new Gaussian state, and the 
detection of the polarization rotation, which is a mea- 
surement of the field observable Xph leads to an update 
of the atomic state, but it retains the Gaussian form [T^ . 
We can hence describe the atomic state by the mean val- 
ues and by the covariance matrix which fully characterize 
a Gaussian state, and we shall now summarize the dy- 
namics of the atoms+fields with an effective formafism 
that describes the time evolution of the mean values and 
the covariances of the atomic state due to interaction and 
measurements. 

We first define a vector of variables (operators) y ~ 
(a^atjPati a;ph,Pph)"^; with the corresponding vector of 
mean values m = (y) and with the covariance matrix 
7 where = 2Re((yj - {yi)){yj - (yj)))- Our Gaussian 
state description of the spin polarized sample and the 
linearly polarized light translates into the specification 
of the initial values, 



m = 

7 = l4x4- 



(4) 
(5) 



and the update formula due to interactions during a time 
interval t 



where 



S = 



m Sm + V 
7 ^ 878^ 

fl 
10 
kV7 1 

\o 1 y 



and 



Let us write 



7 




m = (HiAjms) 



(6) 
(7) 



(8) 

(9) 

(10) 
(11) 



with the covariance matrix for the atomic variables, 
yi = (a^atjPat)^, the covariance matrix for the field 
variables, y2 = {xph,PphV, and the correlation ma- 
trix between yi and y^. When a measurement of the 
variable a;ph is performed, the outcome takes on a ran- 
dom value, given by the Gaussian probability distribu- 
tion. For a short segment of light, the mean value of 
Xph is 777.3 = '«\/Tr(7,2, and the variance is one half (the 
incident field variance is only infinitesimally modified by 
the atoms). The measurement of a value x^eas for a;ph 
collapses the field state and transforms the atomic com- 
ponent according to 



-^7 
rciA 

niB 



A.y — C.Y(7rB-y7r) Cj^ 
niA + C^(7rB^7r)"((a::i, 

l2x2 






m3),0)^ 



(12) 

(13) 
(14) 
(15) 
(16) 



where tt = diag(l,0) and (. .•)^ denotes the Moore- 
Penrose pseudoinverse. To the lowest, relevant order, 
(7rB..,7r)~ = diag(l,0). Equations IH}, and ifTfiji 

"refresh" the atom and field variables corresponding to 
the subsequent light segment which has no correlations 
with the atoms prior to interaction. Note that the quan- 
tity (xmeas — TT^s) IS & random variable with vanishing 
mean and variance 1/2, and since C-y scales with y^r, the 
measurement induced, random displacement of the mean 
values lfT3ll can also be expressed in terms of a Wiener 
increment with zero mean and variance t, cf. Pl, 3. 

Before proceeding to the interaction between the atoms 
and an unknown magnetic field, we note that the dynam- 
ics in is of quantum non demolition type (QND), i.e., 
it permits a detection of the atomic variable Pat without 
changing this observable. Such a detection (carried out 
when the field variable Xph is read out after the atom-light 
interaction) , in effect leads to squeezing of the p^t compo- 
nent of the atomic spin around a random value, which is 
selected by the random outcome of the measurement pro- 
cess and by the deterministic Larmor rotation. Of course 
a single, infinitesimal time step as in Ij2l3|l leads only to 
an infinitesimal squeezing, but as the evolution proceeds 
continuously in time it leads to a monotonic reduction 
of the variance of the atomic spin component. If atomic 
spontaneous decay is taken into account, the squeezing 
has an optimum and further interaction with the optical 
probe leads to incoherent depopulation among the atomic 
states. We are indeed able to describe also such processes 
in the Gaussian state formalism but since the im- 

portance of spontaneous processes is very system specific, 
and it can be reduced by going to sufficiently large detun- 
ings, we shall ignore spontaneous decay in the following. 
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III. ESTIMATION OF AN UNKNOWN TIME 
DEPENDENT B-FIELD 

The i3-field causes the Larmor precession ^ . By prob- 
ing the value of Pat, we acquire information about the 
magnetic field if its value is not already known. In 0,^3, 
we found that a constant magnetic field is effectively 
probed with a time dependent variance on the estimate 
given by 



AB{tf = 



(17) 

with ABq representing our prior knowledge of the field. 
In the limit of n'^t > 1, we have AB{t)'^ ~ 6/(KV^i^) 
which is independent of the prior knowledge and which 
refiects a more rapid reduction with time of the variance 
than expected from a conventional statistical argument. 
This is due to the atomic squeezing, as it progressively 
makes the system more and more sensitive to the mag- 
netic field perturbation. 

In the following, we shall generaHze the analysis to the 
case of time dependent i?-fields. A convenient model for a 
random field is a damped diffusion (Ornstein-Uhlenbeck) 
process, governed by the stochastic differential equation 



dB{t) = --fbB{t)dt + y/^dWb 



(18) 



where the Wiener increment dWb has a Gaussian distri- 
bution with mean zero and variance dt. Our task is to 
expose atoms to a reaHzation of this process and to use 
the polarization measurements to construct an estimate 
for the actual current value of the field. 

The Ornstein-Uhlenbeck process can be simulated on 
a computer, but we can also make statistical predictions, 
e.g., the steady state mean vanishes and the variance is 



Var«t(B) 



2ib 



(19) 



and if the value of the field is estimated at time t^'m the 
laboratory to be Bj^ with a variance Vt^ on the estimate, 
our best estimate for the value at a future time t > II 
takes the value Bt exp(— 7h(t — t^)) with a variance 



F, = 14,e-2^'(*-*-) + -l^(l 
276 



-276(^-tL)^ 



(20) 



The Ornstein-Uhlenbeck process can be very slow, in 
which case the field retains its random value almost con- 
stantly over long times, and we expect to recover the re- 
sults in 01 for estimation of a constant field, because the 
accumulated photo detection record over time carries in- 
formation about the time dependent field, which is known 
to differ not very much. If the Ornstein-Uhlenbeck pro- 
cess is very fast, the accumulated photo detection record 
until the present time t gives only little information about 
the present value. Note that by scaling 7;, and (if, by the 
same factor, the variance ifTfljl is unchanged, but the pro- 
cess changes from slow to rapid fiuctuations. 



Our theoretical description of the estimation process 
deals with a joint Gaussian distribution for the quantum 
variables and the classical magnetic field. As in [a| we for- 
mally treat the B-field as the first component in our vec- 
tor of five Gaussian variables y = (^, a;at,Patj 2:ph,Pph); 
where the tilde is used to distinguish these variables from 
those in the previous section. The Gaussian state is char- 
acterized by its mean value vector in = (y) and its co- 
variance matrix 7 where 7^ = 2Re((yi - (yi))(yj - (yj))). 

Note that although, e.g., x^^t is the same operator in 
this and the previous section, its Gaussian state mean 
value and variance are not the same, because the Gaus- 
sian state is the probability distribution assigned by the 
observer given his or her acquired knowledge about the 
system. In the previous section we determined this prob- 
ability distribution conditioned on full knowledge of the 
time dependent B-field and the photo detection record, 
whereas in the current section, only "the experimental- 
ist's" knowledge of the detection record is assumed. 

By propagating the update formulas for m and 7 
due to the Larmor rotation, the atom-Hght interaction, 
the random outcomes of the probing process, and the 
Ornstein-Uhlenbeck process, we obtain at any time a 
mean value fhi and a variance 711/2 for the i?-field. 

The initial values are 



m = 

7o = diag(2Var(Bo),l,l,l,l) 



(21) 
(22) 



Due to the interaction between the B-field and the atoms 
and between the atoms and the photons we have 



m 



7 



Sim 
S17S 



(23) 
(24) 



where 



Si = 



/ 


1 
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-fiT 





1 


















1 





V 














1 / 



(25) 



In this section the B field is one of the variables, and 
hence the equation Ij23|l is linear unlike the affine trans- 
formation ||SJ. Due to the Ornstein-Uhlenbeck process 
we have 



mi 
7 



mi(l - 
S27S2 



IbT) 

+ L 



(26) 
(27) 



where S2 = diag(l — 7fcT, 0,0,0, 0) and L = 
diag(crbr,0,0,0,0). 

To handle the measurement on Xph we write the co- 
variance matrix and mean value vector as 



7 




m = {m.A,niB) 



(28) 
(29) 
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with A^, the covariance matrix for the i?-field and atoms, 



■7 

(-B,Xat,Pat) 



B™ the covariance matrix for the 



yi 

photons, y2 = 

for yi and y|^. The measurement of Xph then transforms 
these matrices according to 



{xph,Pph) , and the correlation matrix 



— C^(7rB^7r) C'^ 

l2x2 



(ttB^tt) (a;,„ 





(30) 

(31) 

(32) 

(33) 
(34) 



These equations fully describe the conditioned and the 
deterministic evolution of the multi-variable Gaussian 
distribution, and in particular we get access to the es- 
timator for the i3-field in the form of its mean and the 
corresponding covariance matrix element. The input to 
the estimation protocol is the constant parameters of the 
problem and the outcome of the photo detection, which 
will drive the mean values, and hence the estimator, to a 
non-trivial result, cf. Eq. l(33ll . We note that whereas the 
estimator depends on the actual measurement outcome, 
the covariance matrix evolves in an entirely determinis- 
tic manner, and we can hence theoretically predict the 
magnitude of the error as it was also done in Q. 

We wish to simulate the protocol with a given realiza- 
tion of the noisy field, and to this end we have to combine 
the theories of this and the previous section. In such a 
simulation it is the actual field B{t) that acts on the 
atoms and hence leads to the probability distribution for 
the photo detection record, and the measured quantum 
field variable Xph has the mean value ma = Ki/rTO2, given 
by the expressions of the previous section. Hence the 
measurement outcome can be written K^/Tm2+X where x 
is uncorrelated with its value at previous detection times, 
and it has vanishing mean and a variance of 1/2. We are 
thus effectively communicating the value of the simulated 
B-field through the use of the "theorist's" estimator of 
the atomic state in the simulation of measurement out- 
comes. Assuming a random measurement outcome with 
the statistics described, we update the mean value vector 
and covariance matrix of the previous subsection, and we 
use the same simulated measurement outcome in the up- 
date formula l(33jl . The simulated measurement outcome 
may deviate from the value fh^ currently expected by the 
"experimentalist" both because of the noisy contribution 
and because the mean is given by its true value and not 
by his or her estimate. The latter is responsible for driv- 
ing the B-field estimate towards the actual realization in 
our numerical simulations. 



IV. RESULTS 

We now have a complete theory which, given the de- 
tection record, provides the estimator for the time depen- 
dent i?-field and its variance. In practice, the experiment 
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FIG. 2: Uncertainty of the i?-field as a function of 
time. We have used the numerical values of Table 1. 
The dashed/red line is the analytical result for a con- 
stant B-field. The full/black lines represents fluctuating 
fields with the same steady state uncertainty \/ ab/2'yb = 
IpT, but with decreasing rates. Prom above, "ft — 



10'' s-\ 10'*s-\ lO^s" 



,10^ 



For small values of 



7s the full/black curves approach the constant B-field result 
as expected. 



should be run, and the analysis should be applied on the 
full detection record, either simultaneously with the ex- 
periment or afterwards. Some calculation is necessary 
since the -B-field estimator at any time involves knowl- 
edge of the full vector of mean values and the covariance 
matrix. The covariance matrix, and in particular the 
variance of the i3-field estimate evolves deterministically 
with time. We can in fact solve the equations for the 
covariance matrix analytically, and in the long time limit 
we find the steady state variance on our i?-field estimate: 



Var(B) = 



(35) 



When the Ornstein-Uhlenbeck process is slow, early de- 
tection events and estimates provide already an estimate 
for future values of the field, cf. Ij2()|l . which is further 
refined by the continued measurement record. As also 
shown by ll2?1|l . if the rate 76 is high, the estimates quickly 
loose their significance, and only probing for a short time 
before t is useful in the estimate of of B(t). In Fig. [2 
we show how the variance of the i?-field estimate starts 
with the prior steady state value ifTflll before any probing 
takes place, and evolves to the steady state value given 
by II35II ■ In the figure we assume the same steady state 
value for the i?-field variance <7h/2ji, in all curves, but 
with different values of the fluctuation rate constant 76. 
For comparison we show also in the flgure the analytical 
result H17|l for a constant, unknown field with the same 
prior uncertainty. This curve converges to zero, but it is 
clearly seen to follow the curves with finite 7b on short 
time scales (determined by 7^). 

In Fig.|S|we show the result of the application of a given 
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0.1 0.2 0.3 0.4 0.5 -0.2 -0.15 -0.1 -0.05 
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FIG. 3: Time dependent value of the actual, simulated B-field 
(full/red line) and the estimate based on a simulated detection 
record (dashed/green line). The atomic parameters are given 
in Table 1, and the B-field is characterized by 7;, — 10'^ 
and (Ti, = 2 X lO^pT^/s. The inset shows the error of the 
estimate i'Ml as a function of a variable translation of the two 
curves with respect to each other, and the horizontal dashed 
line in the inset shows the error obtained from a weighted 
average of delayed estimates, see text and Eq. I|37|l. 



realization of the Ornstein-Uhlenbeck process. The figure 
shows the actual process as a full line and the estimator 
as a dashed fine for a typical time window. We observe 
that the estimator tracks the gross structure of the de- 
pendence very well, but smaller transients are not repro- 
duced. Taking a second look at the figure, we also ob- 
serve that the estimate is systematically lagging behind 
the true reaHzation of the field. This is not surprising, as 
the estimate makes explicit use of all past measurements 
to predict the current value. In particular, if no measure- 
ment results are provided, the estimate for the value is 
simply obtained from past values according to Eq. itTsjl . 
This should, however, prompt attempts to make an even 
better estimate for the field using not only the detection 
record until the instant of interest, but also the later de- 
tection events. We recall that we are actually probing 
the atomic spin state, and the Larmor precession due to 
the S-field may well be detected also after it took place. 

To justify this increased effort, we have made a very 
simple transformation of the data, consisting in a tem- 
poral displacement of the two curves in Fig. |S[ so that 
we compare the current estimator fhi{t), obtained by our 
theory with earlier values of the field B(t — T). The inset 
of Fig. shows the numerically calculated error of the 
estimate, averaged over a long detection record, plotted 
as a function of the delay T, 

Error2(r) = — [ " {mi{t) - B{t - T)fdt, (36) 

and we indeed see, that for a range of values for T the 
error is smaller than for T = 0, i.e., we obtain a better 
estimate if we assign our estimate to the value of the 
B-field a little earlier. 



FIG. 4: Time dependence of optimal weight factors in 
Eq. (j^ZJ, favoring contributions with a finite delay around 
0.02 ms as suggested by the insert in Fig. |^The atomic pa- 
rameters are given in TableQ] and the B-field is characterized 
by 76 = 10^ s"^ and ab = 2 x 10^ pTVs. 



A slightly more elaborate procedure to improve the es- 
timate obtained from the above procedure is obtained by 
assuming not just a simple delay but a temporal convo- 
lution of the estimators, B^pKi) — a{t — t')rhi{t')dt' . 
Numerically, we have identified the optimum delay dis- 
tribution by minimizing the (squared) error 

^ pto+T 2 

— J [^a^mi{t-idt)- B{t)j dt, (37) 

which turns out to be a linear algebra problem for the 
coefficients a^. Running a series of simulations, we find 
the temporal variation of the delay distribution plotted in 
Fig- El peaking, as expected, around the optimum delay, 
found in Fig. H The error obtained this way is smaller 
than by use of any fixed delay, as indicated by the hori- 
zontal dashed line in the inset in Fig. El 



A. "Gaussian theory of hindsight" 

Both the fixed delay and the weighted average of de- 
lays bring promises for improved sensitivity of magne- 
tometers, if one can wait for the estimate until (a short 
time) after the action of the field. Both procedures, how- 
ever, suffer from their ad hoc character, and in particular 
from the fact that the optimum delay or delay distribu- 
tion are not theoretically available, unless one accepts to 
use simulations as the present ones a guidehne. We shall 
now present the correct theory, which gives the optimum 
estimate and a tight bound on the error without any ad 
hoc procedures. 

The current estimate of the i?-field is given by a Gaus- 
sian probability distribution, and so is the estimate of the 
value at all points in the past, and as measurements on 
the atoms proceed, we may keep improving also our past 
estimate. We hence treat not only the current value but 
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FIG. 5: Time dependent value of the actual, simulated B- 
field (full/red line) and the estimate based on a simulated 
detection record available until 0.1 ms after the action of the 
field (dashed/green line). The atomic parameters are given in 
Table 1, and the B-field is characterized by 75 = 10^ s~i and 
(76 = 2 X 10^ pTVs. 

also past values of the B-field as Gaussian variables to- 
gether with the atom and field variables. To update the 
estimate at an instant T in the past, we need to keep track 
of the entire interval from t — T to the present time t, and 
we extend our Gaussian state formalism by replacing the 
argument yi by a whole vector of values, representing the 
unknown _B-field at discrete times spanning an interval 
from t — T until t. Since we are dealing with the values 
in the past, they are not evolving due to the Ornstein- 
Uhlenbeck process, but inherit their randomly value as 
time proceeds and the elements in the vector of mean 
values and the covariance matrix are simply "pushed to- 
wards the past". Only the current value of B{t) is given 
by the stochastic process. The atomic system is corre- 
lated with both the current and the previous values of 
the field as witnessed by non-vanishing elements in the 
extended covariance matrix, and hence the updating due 
to measurements on the atoms also influence the variance 
of the B-fleld estimate in the past, c.f. the appropriate 
generalization of Eq. ^XM . Now, there is nothing ad hoc 
about the procedure. At any time, we have an estima- 
tor for the value of the fleld over a flnite interval looking 
backwards in time, and we have the variance of these val- 
ues, which is a decreasing function of the time difference, 
approaching a constant, for values so long time ago, that 
their action on the atoms is no longer discernible, and 
hence no further updating takes place due to Eq. 

Fig.0shows a comparison of the actual, simulated field 
with our Gaussian estimator, available 0.1 ms after the 
action of the i?-field on the atoms. Rapid transients are 
not reproduced, but in comparison with Fig. 3, we have 
clearly removed the lag and improved the overall agree- 
ment between the estimator and the actual value of the 
time dependent field. 

Fig.Elreports the variance of the estimate of the B{t)- 
field plotted at a given time in the laboratory after a 
long measurement time, so transients in the atomic dy- 
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FIG. 6: Time dependence of the variance of our esti- 
mate of the time dependent field -B(t)-field around the 
current time tL in the laboratory. The curves represent 
fl uctuatin g flelds with the same steady state uncertainty 
\/ (T6/275 — IpT, but with decreasing rates. Prom above, 
76 = 10*s"\l0^s"i,10^s"Sl0^s"\ The atomic parame- 
ters are given in Table 1. For later times t > tL, we have to 
guess the value from the properties of the Ornstein-Uhlenbeck 
process 12011 . at the current time we have the value 18511 . and 
for earlier times we get the improvement given by the solution 
of the extended Gaussian state updating. The inset shows the 
uncertainty relative to its value at t = tL; improved by an ap- 
proximate factor of 2, when we estimate the past values of 
the field. 



namics have died out. Results are shown for different 
values of but a constant ratio between 7b and ab- The 
estimate for B{t = t^) at the laboratory time is based on 
measurements that have registered the action of the field 
at all previous times and takes the value given by ll35ll . 
For future times t > t^, we have to use the current es- 
timate and extend it by our knowledge of the Ornstein- 
Uhlenbeck process, i.e. by inserting in Eq. lf2(Hl . The 
curves show that it is easier to guess the present value 
than future ones, and that future values are particularly 
hard to guess for a rapidly fiuctuating process. Finally, 
the curves also explore the range oft < tL, where we esti- 
mate past values of the field, and we observe the value of 
hindsight. For small 7^, the uncertainty is small, but we 
see in the inset of the figure, where the variance is plotted 
relative to its value at t — tL, that quite independently 
of the time constant of the fiuctuations, if we can only 
wait long enough, we find an improvement by a factor 
of approximately 2 on the uncertainty compared to the 
equal time estimate. This improvement is a function of 
the physical coupling parameters, and the present study 
suggests a careful analysis of the optimum probing strat- 
egy, including the possibility of time dependent probing 
field strengths and detunings and taking into account also 
atomic decay processes. 
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V. DISCUSSION 

We have presented a formalism, that provides the 
correct estimate of a time dependent B field which is 
known to fluctuate randomly according to an Ornstein- 
Uhlenbeck process. The estimate exhausts the measure- 
ment data and makes the tightest possible conclusions 
from the photo detection record. It is "correct" in the 
sense that any better estimate necessarily requires fur- 
ther information, which could either be in the form of 
prior information about the field or the results of further 
measurements on the system. Our ability at time t to es- 
timate B{t) connects in a natural manner with previous 
results for the estimation of static fields, and it is compa- 
rable with the results of the analysis of time dependent 
fields in Simulations show that this theory actually 
provides a better estimate of the value in the recent past 
(or, instrumentally, if you need to estimate the field at 
time t, you should use the formally obtained estimate at 
time i -I- T for some suitable T) . Our formalism naturally 
generalizes to describe also a scenario, where previous es- 
timates are updated by current measurements, and this 
provides an essential improvement for magnetometers as 
documented in the paper. 

We have not made a complete survey of the optimal 
performance as a function of all physical parameters. At 
this stage, it seems futile to vary all parameters without 
inclusion, in particular, of atomic decay, which will ei- 
ther restrict the magnetometers to finite time analyses 
or which will have to be compensated by a continuous 
optical re-pumping of the atoms. As pointed out in Q, 
the field estimation can be made robust to imprecise in- 
formation about some of the physical parameters, in par- 
ticular the number of atoms which enters the atom-light 
interaction strength. Instead of just calculating the field, 



one can apply, in a feed back set-up, a compensating 
field which should freeze the Larmor precession, inde- 
pendently of the number of atoms, and the field estimate 
is now given by the value of this compensating field. It 
is possible to include a time dependent feed back field in 
our equations, and hence to simulate the performance of 
any feed back strategy and its robustness against fiuc- 
tuations in physical parameters. Another important ex- 
tension of the theory deals with the assumption of the 
Ornstein-Uhlenbeck process with given parameters. The 
parameters enter explicitly in our estimation procedure, 
and if these parameters are not known, or if they are 
functions of time, e.g., as in the case of cardiography on 
a beating heart, more refined theory will bee needed. At 
this point, we recall, that the Ornstein-Uhlenbeck model 
represents our prior knowledge about the field fiuctua- 
tions. A conservative high estimate on the noise fiuctu- 
ations will hence give a similar conservative estimate on 
the field, which might be improved if better limits were 
known. A theoretical investigation of this problem could 
for example simulate the estimation of a noisy field, as- 
suming a different value of the parameters than applied in 
the synthesis of the field, and investigate numerically the 
difference between the variance on the estimate obtained 
from the theory and the actual statistical agreement be- 
tween the estimate and the field. 
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